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Instabilities and turbulence-like dynamics in an oppositely driven 
binary particle mixture 
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Abstract - Using extensive particle-based simulations, we investigate out-of-equilibrium pattern 
dynamics in an oppositely driven binary particle system in two dimensions. A surprisingly rich 
dynamical behavior including lane formation, jamming, oscillation and turbulence-like dynamics 
is found. The ratio of two friction coefficients is a key parameter governing the stability of lane 
formation. When the friction coefficient transverse to the external force direction is sufficiently 
small compared to the longitudinal one, the lane structure becomes unstable to shear-induced 
disturbances, and the system eventually exhibits a dynamical transition into a novel turbulence-like 
phase characterized by random convective flows. We numerically construct an out-of-equilibrium 
phase diagram. Statistical analysis of complex spatio-temporal dynamics of the fully nonlinear 
turbulence-like phase suggests its apparent reminiscence to the swarming dynamics in certain 
active matter systems. 
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Introduction. — Lane formation [T] is one of the rep- 
resentative examples of nonequilibrium phase transitions. 
When two kinds of particles are driven in opposite direc- 
tions, the system exhibits a self-organization from a uni- 
formly mixed state into strongly ordered anisotropic pat- 
terns. This driven segregation phenomenon, first found in 
the computer simulations [2]-[8] , has been observed in labo- 
ratory experiments such as mixtures of oppositely charged 
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or dusty plasmas in presence of an external 
Another important class of examples is 
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pedestrian and traffic flow dynamics 
single-particle motions in colloidal dispersions, a recent ex- 
periment has proposed the underlying mechanism of the 

A lat- 
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lane formation as a dynamical "lock-in" state 
eral mobility of particles is initially enhanced by frequent 
collisions, which however decreases considerably once lane 
is formed. This microscopic dynamics leads to the trap- 
ping of particles within the lanes and thus growth of the 
lane structures. This is also consistent with a physical pic- 
ture employed in the phenomcnological dynamic density- 
functional theory [5j. 

In this letter, using the particle-based simulations, we 
reveal entirely new aspects of such a lane forming sys- 
tem. We show that the anisotropy in the friction force is 
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Fig. 1: (Color online) An illustrative figure of our model. Blue 
particles are driven to the positive x direction while white par- 
ticles are driven to the negative x direction. Each particle 
experiences a friction force proportional to its velocity, whose 
coefficients can be different in x and y directions. 



one of the key parameters that control the onset and sta- 
bility of the lane formation. By systematically changing 
the ratio of two friction coefficients, we find surprisingly 
rich dynamics. More specifically, lanes can become unsta- 
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ble for a sufficiently small transverse/longitudinal friction 
ratio. This transition appears to be triggered by long- 
wavelength interface undulations, which eventually leads 
to turbulence-like convective flows. We numerically con- 
struct a stability diagram as functions of driving speed and 
the friction ratio, and identify the boundary lines sepa- 
rating several qualitatively different dynamical phases. In 
contrast to the lane phase, the structure of turbulence-like 
phase looks considerably isotropic due to random convec- 
tive motions of clusters of the same particle species. How- 
ever, statistical analysis reveals that the spatio-temporal 
order is characterized by correlation lengths that signif- 
icantly exceed the constitutive particle size. Interest- 
ingly, these features are reminiscent of swarming of self- 
propelling elements such as in bacteria colonies 
and flocks of birds or fish 
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Model. — An illustrative figure of the simulation 
setup is shown in fig. [I] Our model system consists of N 
particles interacting with a pair potential <p(r) (explained 
below) on a two-dimensional square plane of length L. The 
particle number density is thus p = N/L 2 . The system is 
driven to out-of-equilibrium state by an externally applied 
field F. To describe how the particles react to F, we as- 
sign the variable s for each particle such that the half of 
the particles with s = +1 are driven to the positive x di- 
rection by F, while the other half with s = — 1 are driven 
to the negative x direction (see fig. [I]). Note that s here is 
introduced as the label only to distinguish the type of re- 
sponses of the particles to F. In contrast to the Brownian 
dynamics studies [2~5|[7], we employ here a deterministic 
particle simulation with inertia, allowing the particles to 
possess arbitrary friction coefficients in the two directions, 
a x and a y , respectively. The friction terms are introduced 
here to mimic the momentum exchange with the substrate. 
The equation of motion for the i-th particle of its position 
Xi(t) for i = 1,2,---, N is 
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where dxi/dt = Vi = (ui,Vi) is the particle velocity, m 
is the particle mass common to all particles. Without 
any interactions, the external force F and the frictional 
force balance in steady state, which drives the particles 
toward their preferential directions at the constant speed 
Vq = F/a x . Hereafter, we use Vq to measure the strength 
of the external driving. To model the short-ranged repul- 
sion between the particles irrespective o/si, we employ the 
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, . < 2 1 / 6 <7 and otherwise </>wca — 0, where r^ — 
We also tested the Lennard- Jones potential 
to evaluate effects of short-ranged attractions, which we 
found are negligible (data not shown). 



Our model given in eq. ([!]) is largely motivated by pedes- 
trian and traffic dynamics in which the self-propelling ele- 
ments interact frictionally with the substrate, rather than 
thermal systems such as colloidal dispersions. The unique 
feature of the present model is the generalization of the 
friction coefficients a x and a y . There are a number of 
examples on the anisotropic frictions such as elongated 
molecules moving in a viscous fluid [23]. In the context 
of pedestrian dynamics, a x and a y are often called "sen- 
sitivity" to the difference between its actual velocity and 
its optimal velocity Vo. Pedestrians may try to keep their 
optimal velocity while they actively change their lanes, 
which implies that a y can be much smaller than a x . 

To minimize effects of the boundaries, we impose the 
periodic boundary conditions for both x and y directions. 
We rescale all variables in the units of energy e, the par- 
ticle diameter er, and the particle mass m, which leads 
to the rescaled time t = tyjma 2 je and the rescaled ve- 
locity V = Vo/^/e/m. We discretize Eq. |l| with the 
time step At, and integrate it by using the velocity Ver- 
let method. To ensure a sufficient numerical accuracy, 
we choose At = 0.001. The initial particle configurations 
are prepared as follows. Particles are first located ran- 
domly on the square lattice points. The magnitude of 
the initial velocity assigned to each particle is fixed to 
be \vi\ = 6y/e/m with random orientations. The initial 
configuration is then obtained by performing 1000 steps 
numerical integrations of Eq. |l]) without introducing any 
dissipation (i.e., a x = a y = 0) in the absence of the exter- 
nal driving (Vq = 0). 



Results. — For the isotropic friction case a y /a x = 1, 
the uniformly mixed initial state becomes unstable and 
the usual lane formation proceeds. In contrast, for small 
transverse friction case a y /a x -C 1, this lane formation is 
only transient and is unstable against long-wavelength un- 
dulations of the interfaces. In fig. [2] typical snapshots and 
the corresponding velocity field obtained from the simula- 
tions for a y /a x = 0.2 are shown. A transiently formed lane 
(fig. [2] (a)) becomes unstable and buckle (fig. [2] (b)), lead- 
ing to the turbulence-like convective flow structure (fig. [2] 
(c)). The velocity fields shown on the right hand panel of 
fig.[2](c) appear to be very complicated, where the complex 
transient jet or vortex structures are observed. Statisti- 
cal properties of this turbulence-like phase will be further 
investigated later on. To our knowledge, this anisotropy- 
driven instability was previously unknown, and is the first 
main result of our present study. While we performed the 
simulations for increasing system sizes, L/a = 40,80 and 
120 with the fixed particle number density (per 2 = 1.0), no 
qualitative and systematic changes were observed. We will 
therefore show the results obtained from the simulations 
for L = 80a and N = 6400 (i.e., pa 2 = 1.0) below. 

To characterize the observed patterns more quantita- 
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Fig. 2: (Color online) Typical time evolution of a transiently 
formed lane structure in its unstable regime observed in our 
simulations. The transition to a turbulence-like flow proceeds 
from (a) to (c). The parameters used in (a)-(c) are Vo = 5.0, 
a y /a x — 0.2, L = 80cr and pa 2 = f.0. Left panel: snapshots 
of the particle configurations. Right panel: the corresponding 
velocity field. Arrows show the orientation of the velocity and 
colors indicate the magnitude of the velocity. 



tively, we introduce the (local) temperature Ta as 
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where O is the time average of O between t = 5000 
and 10000 for any variable O. A denotes a horizontal 
stripe area of a x L in the simulation box and (O,) A — 
12ieA Q f° r anv variable O,-. Here Na is the number 
of particles (white+blue) inside A. 

Note that Ta measures the kinetic energy in the local 
co-moving frame with its local mean velocity (v) , and is 
higher in the regions where collisions are more frequent 
such as in the interfaces between oppositely moving lanes. 
In fig. [3] (b) , we show the spatial profile along y direction 
of the local particle number density pa — N^/(aL), the 
temperature Ta and the x component of the mean velocity 
(it) as a function of y in the stable laning phase obtained 
for a y /a x = 1.0, V = 4.0, L = 80er and pa 2 = 1.0. 




Fig. 3: (Color online) A typical snapshot (left panel of (a)) 
and the corresponding velocity field (right panel of (a)) of lane 
configuration observed in our simulations. Profiles of the total 
(white+blue) particle number density pa (left panel of (b)), the 
temperature Ta (also left panel of (b)) and the x component 
of the mean velocity (right panel of (b)) in the stable lane 
forming phase as functions of y. Simulation parameters used 
are a y /a x — 1.0, Vb = 4.0, L = 80ct and pa 2 — 1.0. 

The temperature Ta of the lane phase (shown in fig. [3] 
(b)) is almost zero inside the lanes, whereas it is higher 
in the boundary layers between the oppositely oriented 
lanes. In contrast, the number density pA is lower in the 
boundary layers and is higher inside the lanes. This is be- 
cause there are much more collisions between the particles 
and considerably disordered particle configurations are re- 
alized within the layers. For a sufficiently small a y /a x , the 
collision-induced shear stress at the boundary layer would 
be enough to induce collective migrations of the particles 
transverse to the flow direction, which ultimately causes 
the buckling 24 ( 25 of the lanes (see fig [2]) . 



Phase diagram. — In addition to the lane and 
turbulence-like phases, we also find a few other character- 
istic phases depending on Vb and a y /a x , such as the com- 
pletely jammed state in which clusters of particles push 
against each other without any motion nor transport of 
particles. To quantify the various patterns observed in 
our systems, we introduce the order parameter $ as fol- 
lows. First, we slice our simulation box into thin sections 
parallel to x axis of width a. Let nf and n~ be the number 
of particles with s = +1 and s = — 1, respectively, found 
in the i-th section, then we define our order parameter 



N s ^ nj 



(4) 



where N s is the total number of the thin sections, and the 
bracket represents the long time average in steady state. 
By definition, $ changes from to 1. For randomly mixed 
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Fig. 4: (Color online) The out-of-equilibrium phase diagram 
that is numerically constructed by looking at the order param- 
eter $ defined in eq. Q, on the plane spanned by the external 
driving speed Vb and the friction constant ratio a y /a x . The 
number density and the system size are fixed as pa 2 = 1.0 and 
L — 80a. The symbols on the diagram indicate the different 
values of $: O: $ is from 0.00 to 0.25; ♦: <E> is from 0.25 to 
0.50; □: $ is from 0.50 to 0.75; ■: $ is from 0.75 to 1.00. 
Symbol + and X indicate non-flow phases: +: randomly mixed 
initial configuration; X: completely jammed state. Note that 
the value of $ in the non-flow phase is around 0.1. 



configurations, $ is close to zero, while it approaches one 
as the system evolves the laning order. We note that our 
definition of 4> is essentially the same as those introduced 
previously [2 ., and small differences are insignificant for 
describing the qualitative aspects of our results. 

With the order parameter $ defined above, we now con- 
struct a phase diagram. In fig. |3J the values of $ ob- 
tained from our simulations are shown on the Vq — (a y /a x ) 
plane. Five qualitatively different phases, the randomly 
mixed, stable lane, jammed, uniform gas, and convective 
or turbulence-like phases, are identified for a fixed num- 
ber density pa 2 — 1.0. When Vq is smaller than 2.0, the 
external force is insufficient to produce any macroscale 
particle flow, where the system almost keeps its initial 
configuration. Once Vq becomes larger than 2.0, particles 
start to move and the large-scale pattern formation pro- 
ceeds. For Vq slightly larger than its critical value 2.0, 
the usual lane forming dynamics is observed if the dissi- 
pation is almost isotropic, for which we obtain $ rs 1. As 
Vq further increases, the kinetic temperature T increases 
significantly, leading to a uniform gas-like phase. This 
process is characterized by a gradual broadening of the 
lane interfaces, and we note that this gas-like phase was 
addressed previously in Ref. j6j. However, the lane desta- 
bilization observed in our system is distinct from this gas 
phase, because it is induced only by the strong anisotropy 
in the friction coefficients. Also, the resulting turbulence- 
like phase is characterized by the collective movements 
of the particles, rather than the gas-like mixture. As Vq 




Fig. 5: (Color online) Time evolution of the order parameter 
$ during the extensive long time simulation for Vb = 5 and 
a y /a x — 0.5. In this case, intermittent lane formations are ob- 
served between laning phase and the turbulence-like convective 
flow phase. 

increases, $ decreases monotonically for a y /a x 1, but 
for smaller a y /a x (typically less than 0.5), $ changes al- 
most discontinuously from the lane to the turbulence-like, 
from the lane to the jammed, and from the jammed to the 
turbulence-like phases (see Fig. 4). While the former is the 
continuous crossover, the latter may arise from a certain 
hydrodynamic-type instability. 

At several points indicated by "X" in fig. [4j $ is excep- 
tionally small as compared with those of the neighboring 
points. These points correspond to the jammed phase, 
where the particles do segregate but their interface ap- 
pears perpendicular to the driving direction. In this case, 
the whole system maintains its static force balance state 
without any particle currents. As Vq is increased, the ver- 
tical interface becomes unstable against long wavelength 
undulations, and the particles start to move. 

Turbulence-like phase- As a y /a x is decreased, the sys- 
tem begins to show complex spatio-temporal dynamics. 
For example, at a point indicated by * (as well as its small 
surrounding region) on the diagram in fig. |4j we observe 
a unique "quasi-periodic" time evolution of $, which is 
shown in fig. [5] The lane structure (with $«1) breaks up 
stochastically, but the resulting disorganized state ($ « 0) 
is unstable and the lane formation immediately proceeds 
again, resulting in the spikes of 4> in fig. [5j Strikingly, al- 
though the precise mechanism of this nonlinear oscillation 
is at present unclear, this whole process was observed to 
repeat during our extensive long time simulations, and not 
transient. 

For even smaller values of a y /a x , the system is brought 
into a highly chaotic or turbulence-like flow phase, which 
contrasts to the highly ordered and anisotropic lane phase. 
To quantitatively characterize this phase, we look at the 
correlation function of y-component instantaneous veloc- 
ity field, which is defined as 



C yy (x, y) 



(v( x o> Vo? t)v(x + x , Vo + V, t)) 
(v(x ,yo,t)v(x ,yo,t)} 



(5) 



where (• • •) represents average over time t and space Xq, yQ. 
Whereas the external force is absent in y direction, the cor- 
relation function is strongly anisotropic, as seen from the 
inset of fig. [6] This feature is also understood by compar- 
ing the profiles of C yy along the x-axis (with y = 0) and 
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the time derivative of the total energy E = dE/dt as 
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Fig. 6: (Color online) Spatial correlation function of the y- 
component of the instantaneous velocity field obtained from 
the simulations for a y /a x — 0.0 and Vo = 4.0. In the main 
panel, C yy (x,y = 0) (red and solid line) and C yy (x = 0,y) 
(green and dashed line) are plotted. Inset figure shows the 
two-dimensional density plot of C yy to highlight the anisotropic 
feature of the velocity correlations. 



along the y-axis (with x = 0), shown in the main panel of 
fig. [6] The transverse correlation C yy (x,y = 0) becomes 
negative at around x c /a f=a 16, while the longitudinal one 
C yy (x = 0, y) decays exponentially from the origin and be- 
comes slightly negative at around y c / 'a ss 30. This obser- 
vation provides us a physical picture that densely packed 
particles move in a coordinated manner and form an elon- 
gated cluster or jet of typical sizes x c and y c . In addition, 
the clusters of the oppositely moving directions encounter 
in several places, as the upward and downward flow are 
neighboring in fig. [2j which gives rise to transient vortex 
structures. While external driving is solely along the x 
direction, the particles find their escape ways to the less- 
dissipating y-direction. Thus the energy is transmitted 
via the anisotropic dissipation from the x to y direction, 
which gives rise to the distinct collective dynamics at in- 
termediate length scales (see fig. [6]) . A rich variety of the 
extended spatio-temporal dynamics in our lane forming 
system is the second main result of this paper. 



Discussion. — Here we outline simple energetic ar- 
guments for the stability of lane structures, which might 
provide some (but rough) physical interpretations on the 
observed instabilities. The total energy of our system con- 
sists of the kinetic and potential energies and is written 
as 



N 



mvt 



+ 
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(6) 



Now we introduce the drift velocity — SjUj 
and fluctuations from it as Sui = Ui — s^Ud- Note that 
8vi — Vi because we can obtain from eq. |l| that YliLi v i = 
in steady state. With these quantities, we can compute 



E = Na x Ud(V Q - ltd) 



N 



a x (5ui) 2 + a y (Svi) 2 ] . (7) 



The first term on the right hand side of eq. |7|, W = 
Na x VoUd, represents the work done by the external force 
per unit time. The second term, Qhk = ^dxu\, known 



as the rate of housekeeping heat 26 , corresponds to the 



average energy dissipation necessary to keep the whole 
system at this steady state. Third terms in eq. 
Qex — Hj = i [ a x{3ui) 2 + a y (5vi) 2 ] , represents the "ex- 
cess" energy dissipation per unit time, which arises from 
fluctuations of particle velocities from its steady-state av- 
erage. Equation Q suggests the usual energy balance (per 
unit time), E = W — Q, where Q — Qhk + <3cx- 

Whereas our model does not have any direct route 
to thermal equilibrium, the particle motions are under- 
damped and the system is expected to be fully thermalized 
at a driven steady state owning to their inertia. Given the 
principle of minimum entropy production in such driven 
steady state (or something analogous to this) , the particle 
configurations that can minimize Q cx should appear as the 
most stable one. When both a x and a y have large enough 
positive values, it is achieved by minimizing 5ut and 5ui, 
which indicates that the system prefers a stable lane for- 
mation. However, when a y becomes sufficiently small, ve- 
locity fluctuations along y direction can be unboundedly 
large, which could ultimately lead to a buckling of lane 
structures. Whereas the validity of the minimum entropy 
production principle to our system is entirely unclear at 
present Q such a physical picture is fully consistent with 
what we have found in the simulations and could be useful 
to formulate a more elaborate theory in future studies. 

Finally, we briefly remark on the turbulence-like phase. 
This dynamics is governed by the highly coordinated mo- 
tions of the clusters of the same particle species. The 
instantaneous velocity field exhibits complex vortex and 
jet structures, whose typical sizes considerably exceed 
the constitutive particle size (see fig. [2]). These features 
present some similarities with those of swarming phase 
found in bacterial suspensions or colonies 18 19 . In ac- 



tive matter systems, however, there are uptake of free en- 
ergy and the systematic movement under the "force-free" 
condition at the level of each active elements [16] . Such a 
"micro-to-macro" energy transport is characteristic in ac- 
tive matter, and obviously contrasts to our lane forming 
system where the energy input is solely at the macroscopic 
scale. Nevertheless, this apparent similarity might suggest 
a certain underlying mechanics that is common to driven 
soft matter and active matter systems. Elucidating the 
microscopic particle dynamics in our system would pro- 
vide more quantitative insights into the observed instabil- 



1 We note that there is no universal principle for the pattern selec- 
tions in dissipative structures even when a system is close to (ther- 
mal) equilibrium. See for example ref. [271. 
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ity, which will also help develop the linear and nonlinear 
stability analyses. 

Summary. — Our study opens up the entirely new as- 
pects of the lane forming system, i.e., the prototype model 
for non-equilibrium phase transitions, with a wealth of 
novel complex dynamical behavior. The results obtained 
here will be potentially important not only in technologi- 
cal applications such as microfluidics, but also for under- 
standing collective human behavior such as crowd controls 



or panic dynamics in densely populated areas 15 
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